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Abstract 

WRKY proteins comprise a large family of transcription factors that play important roles in plant defence regulatory 
networks, including responses to various biotic and abiotic stresses. To date, no large-scale study of WRKY genes 
has been undertaken in grape (Vitis vinifera L.). In this study, a total of 59 putative grape WRKY genes [VvWRKY) were 
identified and renamed on the basis of their respective chromosome distribution. A multiple sequence alignment 
analysis using all predicted grape WRKY genes coding sequences, together with those from Arabidopsis thaliana and 
tomato [Solanum lycopersicum), indicated that the 59 VvWRKY genes can be classified into three main groups (l-lll). 
An evaluation of the duplication events suggested that several WRKY genes arose before the divergence of the grape 
and Arabidopsis lineages. Moreover, expression profiles derived from semiquantitative PGR and real-time quantitative 
PGR analyses showed distinct expression patterns in various tissues and in response to different treatments. Four 
VvWRKY genes showed a significantly higher expression in roots or leaves, 55 responded to varying degrees to at 
least one abiotic stress treatment, and the expression of 38 were altered following powdery mildew (Erysiphe necator) 
infection. Most VvWRKY genes were downregulated in response to abscisic acid or salicylic acid treatments, while 
the expression of a subset was upregulated by methyl jasmonate or ethylene treatments. 

Keywords: Evolution, expression profile analysis, grape {Vitis vinifera L.), phylogenetic analysis, synteny analysis, WRKY genes. 



Introduction 

Transcription factors are proteins that bind to specific 
DNA sequences in the promoter regions of genes, thereby 
regulating their transcription. Consequently, transcription 
factors play pivotal roles in numerous plant signalling and 
regulatory networks (Hwang et al, 2010). WRKY proteins, 
which are characterized by a highly conserved domain of 
about 60 amino acid residues, comprise a class of transcrip- 
tion factors that are known to function as transcriptional 
activators or repressors in a number of developmental and 
physiological processes (Eulgem et al. , 2000; Rushton et al. , 



2010, 2012). The first WRKY gQ^Q to be cloned and char- 
acterized was from sweet potato (Ishiguro and Nakamura, 
1994) and this was followed by studies of WRKY gtat?, 
from Arabidopsis (Eulgem et al, 2000), rice (Wu et al, 
2005), and barley (Mangelsen et al., 2008). Moreover, 
large-scale genome- wide studies of WRKY gtaQ?, have been 
described for Arabidopsis (Dong et al, 2003), rice (Ryu 
et al, 2006), poplar (He et al, 2012), tomato (Huang et al, 
2012), cucumber (Ling et al, 2011), and coffee (Ramiro 
et al, 2010). 
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The two most prominent and defining structural charac- 
teristics of WRKY proteins are the WRKY domains and 
the zinc-finger motifs (C-X4_5-X22_23-H-Xi-H or C-X7-C- 
X23-H-X1-C), which provide the basis of their classification 
into groups I-III (Eulgem et al., 2000). Proteins from group 
I contain two WRKY domains and a C2H2 zinc-finger motif, 
while those from groups II and III only contain one WRKY 
domain and a C2H2 or C2HC zinc-finger motif, respectively 
(Eulgem et a!., 2000; Zhang and Wang, 2005). WRKY pro- 
teins bind to the W-box ((C/T)TGAC(T/C)) of their tar- 
get genes (Ciolkowski et al., 2008) as well as a cw-element 
(SURE), which plays a role in promoting transcription (Sun 
et al., 2003). 

WRKY transcription factors are involved in regulatory 
processes mediated by various biotic and abiotic stresses 
(Eulgem and Somssich, 2007; Pandey and Somssich, 2009). 
For example, defence responses in Arahidopsis to the bacte- 
rial pathogen Pseudomonas syringae have been associated 
with the action of AtWRKY38 and AtWRKY62 (Journot- 
Catalino et al. , 2006; Kim et al. , 2008) and the involvement 
of WRKY gQTLQ?, in the resistance of grapevine to pathogens 
was suggested by studies using transgenic tobacco fines 
expressing VvWRKYl, reflecting the important function of 
WRKY gQXiQi, in biotic stress tolerance (Mzid et al, 2007). 
WRKY gQQQ?, are also related to abiotic stress induced gene 
expression: AtWRKY25 and AtWRKY33 show altered 
expression following either heat (Li et al., 2009; Li et al., 
20 11) or salt treatments (Jiang and Deyholos, 2009) and the 
expression of TcWRKY53 from alpine pennygrass (Thlaspi 
caerulescens) is affected by cold, salt, and polyethylene 
glycol treatments (Rushton et al, 2010). Studies involving 
other plants, including rice (Shimono et al, 2011), wheat 
(Wu et al, 2008) and poplar (He et al, 2012), have simi- 
larly provided insights into the diversity of WRKY gene 
function. 

Grape {Vitis vinifera L.) is used both as a fresh food 
commodity and for processed food products such as juice, 
raisins, and wine. It is cultivated worldwide and has great 
economic value (Kreamer, 1995; Jaillon et al, 2007) and so 
there is considerable interest in identifying genes to improve 
grape horticultural characteristics, such as those that pro- 
mote stress resistance. In this regard, the various studies of 
WRKY gQnQi, from different plant species mentioned above 
suggest that members of the grape WRKY gene family 
{VvWRKY) have considerable potential for contributing to 
aspects of stress resistance. This work reports the identifica- 
tion of 59 putative Vv WRKY genes, together with an analy- 
sis of their exon-intron organization and associated gene 
duplication events in the context of gene evolution, since it 
has been shown that gene dupHcation has played a critical 
role in Arahidopsis and rice WRKY gene family expansion 
(Cannon et al, 2004). In addition, this work determined the 
expression profiles of Vv WRKY genes, in six different tis- 
sues and measured their transcript abundance in response 
to different phytohormone treatments and under various 
abiotic and biotic stresses. This study provides a foundation 
for future studies of VvWRKY gene family evolution and 
function. 



Materials and methods 

Identification and annotation of grape WRKY genes 

An Hidden Markov Model profile of the WRKY DNA-binding 
domain (PF03106) was downloaded from the Pfam protein family 
database (http://pfam.sanger.ac.uk/) (Finn et al, 2010) and used 
to identify putative WRKY genes/proteins from the grape genome 
sequence (http://www.genoscope.cns.f) (Jaillon et al, 2007) using the 
BLASTP program and default parameters (Ling et al, 2011). All 
non-redundant gene sequences encoding complete WRKY domains 
were selected as putative WRKY genes. Expressed sequence tags 
(ESTs) of each gene sequence from the public V. vinifera EST data- 
base were used for further validation. The identified grape WRKY 
genes were annotated based on their respective chromosome distri- 
bution (Ling era/., 2011). 



Multiple sequence alignment, phylogenetic analysis, and 
classification of grape WRKY genes 

A total of 59 predicted VvWRKY proteins, with amino acids span- 
ning the WRKY core domain, were included in multiple sequence 
alignments using CLUSTALX version 2.0. 12 (Larkin et al , 2007) and 
Boxshade (http://www.ch.embnet.org/software/BOX_form.html). 
A further multiple sequence alignment including VvWRKY genes and 
those from Arahidopsis {At WRKY) and tomato (Solarium lycopersi- 
cum, SIWRKY) was performed using CLUSTALW. A phylogenetic 
tree based on the alignment was constructed using MEGA 5.0 with 
the neighbour-joining method and with the bootstrap test replicated 
1000 times (Tamura e?a/. , 201 1). Based on the multiple sequence align- 
ment and the previously reported classification of At WRKY genes, 
the VvWRKY genes were assigned to different groups and subgroups. 

Exon-intron structure, tandem duplication, and synteny analysis 
of grape WRKY genes 

The exon-intron structures of the grape WRKY genes were deter- 
mined based on alignments of their coding sequences and their 
respective full-length sequences (Grape Genome Browser; http:// 
www.genoscope.cns.fr/externe/GenomeBrowser/Vitis/), while dia- 
grams were obtained from the online program Gene Structure 
Display Server (GSDS: http://gsds.cbi.pku.edu.ch). WRKY genes 
with tandem duplication events were defined as adjacent homolo- 
gous genes on a single chromosome, while gene duplication events 
between different chromosomes were characterized as segmental 
duplications (Liu et al. , 201 1 ). The specific physical location of each 
VvWRKY gene on its individual chromosome therefore determined 
whether it was regarded as a genes resulting from a tandem duplica- 
tion event. The syntenic blocks used for constructing a synteny anal- 
ysis map of the grape WRKY genes, as well as between grape and 
Arahidopsis WRKY genes, were obtained from the Plant Genome 
Duplication Database (Tang et al , 2008) and the diagrams were gen- 
erated by the program Circos version 0.63 (http://circos.ca/). 



Plant material and treatments 

Grape organs/tissues (roots, tendrils, leaves, inflorescences, fruit, 
and stems) were obtained from 2-year-old 'Kyoho' ( Vitis labrusca 
X V. vinifera) grape seedlings, which were grown in 12 dm^ pots in 
the greenhouse of Northwest A&F University, Yangling, Shaanxi, 
China (34° 20' N 108° 24' E). The third to fifth fully expanded young 
grapevine leaves beneath the apex were taken for hormone treat- 
ments, at which time the shoots of the vines were 25-35 cm long. 

Hormone treatments were carried out by spraying leaves with 
300 nM abscisic acid (ABA), 100 ^iM salicylic acid (SA), 50 ^iM 
methyl jasmonate (MeJA), and 0.5 g/1 ethylene (Eth), and then leaves 
were sampled at 0.5, 1, 6, 24, and 48 h post treatment (Li et al, 
2010). Leaves sprayed with sterile water and similarly harvested were 
used as a negative control. A sahnity stress treatment was carried 
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out by irrigating plants with 2 1 of 250 mM NaCl in the pots (Upreti 
and Murti, 2010) followed by sampling leaves at 1, 6, 24, and 48 h 
post treatment. Seedlings irrigated with 2 1 tap water were used as a 
negative control. A drought stress treatment was performed by with- 
holding water (Yang et a!., 2012) followed by samphng at 24, 48, 96, 
144, and 168h post treatment. Plants were rewatered after 168h of 
drought stress and sampled again 48 h later. Grape seedlings grown 
without drought stress were used as a control. 

Samples of the powdery mildew fungus {Erysiphe necator) were 
used to inoculate young leaves of V. quinquangularis 'Sliang-24' 
(Northwest A&F University, Yangling, Shanxi, China; Wang et al. , 
1995). Leaves were harvested at 6, 12, 24, 48, 72, 96, and 120h post 
inoculation and uninoculated leaves served as a negative control. 

At each time point of each treatment, six leaves from six separate 
plants were combined to form one sample, and all of the treatment 
experiments were performed in triplicate. All these plant samples 
were immediately frozen in liquid nitrogen and stored at -80°C until 
RNA extraction. 



Semiquantitative reverse-transcription PCR anaiysis 

Total RNA samples were extracted from leaves using the EZNA 
Plant RNA Kit (R6827-01, Omega Bio-tek, USA). First-strand 
cDNA was synthesized by reverse transcription of 500 ng total RNA 
using PrimeScript RTase (TaKaRa Biotechnology, Dalian, China). 
The concentration of the cDNA was adjusted using PCR and the 
grape Actinl gene (GenBank accession number AY680701) with 
the primers F (5'-GATTCTGGTGATGGTGTGAGT-3') and R 
(5'-GACAATTTCCCGTTCAGCAGT-3'). Gene-specific primers 
for each VvWRKY gene were designed using Primer Premier 5.0 and 
optimized using oligo 7 (Supplementary Table SI available at JXB 
onhne). Semiquantitative reverse-transcription (RT) PCR reactions 
were conducted using the following profile: initial denaturation at 
94°C for 2min, followed by 30^0 cycles of denaturation at 92°C for 
30 s, annealing at 60 ± 5°C for 30 s, extension at 72°C for 30 s, and final 
extension at 72°C for 2min. PCR products were separated on a 1.5% 
(w/v) agarose gel with ethidium bromide staining and imaged under 
UV light for further gene expression analysis. Each reaction was 
repeated three times and the three independent analyses showed the 
same trends for each gene and treatment. The expression data from 
the semiquantitative RT-PCR were collated, analysed, and visualized 
using the programs GeneSnap and Mev 4.8.1 (Saeed et al, 2006). 

Reai-time quantitative PCR 

Real-time quantitative PCR was conducted using SYBR green 
(TaKaRa Biotechnology) on an IQ5 real time-PCR machine (Bio- 
Rad, Hercules, CA, USA) with a final volume of 20 |il per reac- 
tion. Each reaction mixture contained 10.0 |il SYBR Premix Ex 
Taq II (TaKaRa Biotechnology), 1.0 |.i,l cDNA template, 0.8 |xl each 
primer (1.0 ^iM), and 7.4 |j.l sterile distilled II2O. Each reaction was 
performed in triplicate. Cycling parameters were 95°C for 30 s, 40 
cycles at 95°C for 5 s, and 60°C for 30 s. Melt-curve analyses were 
performed using a program with 95°C for 15 s and then a constant 
increase from 60°C to 95°C. Gene-specific DNA primers were the 
same as those used for semiquantitative RT-PCR (Supplementary 
Table SI). The grape Actinl gene was used as the internal reference 
gene. The software IQ5 was used to analyse the relative expression 
levels using the normalized-expression method (Hou et al., 2013). 



Results 

Identification of WRKY genes in tlie grape fV. vinifera 
L) genome 

A total of 61 genes were originally obtained between 
PF03106 (a Hidden Markov Model profile of WRKY 



DNA-binding domain) and Grape Genome (12X) with 
BLASTP. Based on the presence of apparently complete 
WRKY domains, 59 genes were subsequently selected 
and annotated as being grape WRKY genes. Genes with- 
out a complete predicted WRKY domain were removed 
(GSVIVT01016690001, GSVIVTO 103 140 1001). Further 
analysis of the protein sequences using the NCBI website 
revealed that GSVIVT01016690001 shares 63% similarity 
with a Theobroma cacao GDSL-motif lipase protein, indi- 
cating that GSVEVTO 10 16690001 may not belong to the 
WRKY family. The other putative 59 grape WRKY genes 
were mapped onto the 19 grape chromosomes and then 
renamed from VvWRKYl to VvWRKY59 based on their 
distributions and relative linear orders among the respec- 
tive chromosome. VvWRKY4 (GSVIVT01001332001) was 
located to chromosome l_random region and VvWRKY59 
(GSVIVTO 1007006001) in an unknown region, and were thus 
unlike the other 57 VvWRKY gQwe?,. At least one EST for 52 
VvWRKY gQnQS was found in the pubHc V. vinifera EST data- 
base at NCBI and used for further corroboration, while only 
seven (VvWRKY9, VvWRKYl 3, VvWRKY 14, VvWRKY 
17, VvWRKY 24, VvWRKY 41, VvWRKY 51) were without 
a corresponding EST. However, in order to perform a sys- 
tematic and comprehensive analysis of all of the VvWRKY 
genes, 59 specific primers were designed (Supplementary 
Table SI) for the expression analyses and gene confirmation. 
Detailed information about each VvWRKY gme is showed in 
Table 1, including the WRKYgsm group numbers, gene locus 
numbers, accession numbers for the full-length sequences at 
NCBI, chromosome distribution (start sites and end sites), 
and the length of coding sequences. 

Mu/f/p/e sequence alignment of VvWRKY genes 

The most prominent structural feature of WRKY proteins is 
the WRKY domain, which has been shown to interact with 
the W-box (C/T)TGAC(T/C), thereby activating a large num- 
ber of defence-related genes (Eulgem et al. , 2000). The WRKY 
domain consists of a highly conserved hepta-peptide stretch 
of WRKYGQK at the N-terminus, followed by a zinc-finger 
motif (Eulgem et al. , 2000). A multiple sequence alignment of 
the core WRKY domain, spanning approximately 60 amino 
acids of all 59 VvWRKY proteins is shown in Fig. 1. A total 
of 54 VvWRKY pwtdm were found to have the highly con- 
served sequence WRKYGQK, while the others (VvWRKYS, 
VvWRKYl 3, VvWRKYM, VvWRKYl?, Vv WRKY2 4) vary hy 
a single amino acid. Of these WRKYGKK is the most com- 
mon, consistent with studies in tomato, being present in four of 
the five variants, while VvWRKYl? contains a WKKYGQK 
sequence. One possible consequence of variation in this 
WRKY domain is an altered binding specificity in the DNA 
targets, but this remains to be demonstrated. As previously 
described (Eulgem et al. , 2000), the metal-chelating zinc-finger 
motif (C-X^5-X22-23-H-Xi-H or C-X7-C-X23-H-X1-C) is 
another important stiuctural characteristic of WRKY proteins. 
However, some incomplete zinc-finger motifs were also identi- 
fied, including examples encoded by VvWRKY?, VvWRKYl?, 
and VvWRKY46. Interestingly, in contrast to group-Ill WRKY 
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Table 1 . Grape WRKY genes and accessions 

CDS, coding sequence; NG, no group; NSG, no subgroup, ORF, open reading frame. 
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VvWRKY59 
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29 699 639 
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VvWRKY46 


GSVIVT01 Oil 472001 


CBI22264.3 
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889 


14 


29 916 400 


29 925 511 


9112 
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GSVIVT01 01 4854001 


CBI39865.3 
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1 0 665 036 


1 0 669 055 
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16 666 741 
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CBI36524.3 


1500 


499 


11 


7 835 815 


7 846 1 37 
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VvWRKYIS 


GSVIVT01 024624001 


CBI1 5865.3 


1713 


570 
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8 290 025 


8 295 109 


5085 
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GSVIVT01 025562001 


CBI32633.3 


1317 
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1 4 039 562 


6112 
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CBI28412.3 


1095 
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12 


9 116 731 


9 122 807 


6077 
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GSVIVTOI 030258001 


CBI1 8092.3 


1542 


513 


8 


9 796 056 


9 798 907 


2852 
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GSVIVT01 035965001 


CBI21 139.3 


1593 


530 


4 


6 569 931 


6 576 637 
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19 
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16 094 219 


1 6 096 253 
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605 
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6316 168 
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GSVIVTOI 020060001 


CBI32009.3 


1785 


594 
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10 977 206 


10 982 423 
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VvWRKY22 


GSVIVTOI 028244001 


CBI37053.3 


1440 


479 


7 


4 899 927 


4 903 1 75 


3249 
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VvWRKY40 


GSVIVTOI 029688001 


CBI1 9998.3 


1473 


490 


12 


13 065 135 


13 099 628 
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CBI28821.3 


1497 


498 


12 


5 678 707 


5 681 171 


2465 
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CBI26664.3 


1491 


496 


19 


6 882 419 


6 884 988 


2570 


0 


VvWRKY53 


GSVIVTOI 008553001 


CBI15677.3 


456 


151 


17 


922 600 


925 171 


2572 


0 


VvWRKY3 


GSVIVTOI 01 0525001 


CBI27681.3 


570 


189 


1 


21 460 123 


21 461 397 


1275 


c 


VvWRKYI 


GSVIVTOI 01 21 96001 


CBI27268.3 


852 


283 


1 


628 682 


633 595 


4914 


0 


VvWRKY47 


GSVIVTOI 01 8300001 


CBI1 6682.3 


687 


228 


15 


1 1 498 970 


1 1 504 729 


5760 


c 


VvWRKY37 


GSVIVTOI 020864001 


CBI22108.3 


936 


311 


12 


878 934 


881 289 


2356 


0 


VvWRKY33 


GSVIVTOI 021 397001 


CBI30827.3 


960 


319 


10 


4 894 476 


4 896 340 


1865 


c 


VvWRKY24 


GSVIVTOI 022245001 


CBI21 522.3 


582 


193 


7 


17 794 379 


17 797 240 


2862 


c 


VvWRKY25 


GSVIVTOI 022259001 


CBI21 534.3 


681 


226 


7 


17 958 306 


17 960 930 


2625 


c 


VvWRKY49 


GSVIVTOI 026969001 


CBI40411.3 


606 


201 


15 


18 940 954 


18 942 146 


1193 


c 


VvWRKY21 


GSVIVTOI 0281 47001 


CBI36970.3 


909 


302 


7 


4 200 1 60 


4 202 241 


2082 


c 


VvWRKY44 


GSVIVTOI 033063001 


CBI21329.3 


549 


182 


14 


25 479 103 


25 481 683 


2581 


c 


VvWRKY13 


GSVIVTOI 0331 94001 


CBI24019.3 


471 


156 


4 


9 399 944 


9 400 803 


860 


0 


VvWRKYU 


GSVIVTOI 0331 95001 


CBI24020.3 


306 


101 


4 


9 409 805 


9 41 1 286 


1482 


c 


VvWRKYW 


GSVIVTOI 034968001 


CBI22862.3 


930 


309 


5 


530 084 


531 773 


1690 


c 


VvWRKYS 


GSVIVTOI 035426001 


CBI20684.3 


501 


166 


4 


1 209 585 


1 211 712 


2128 


d 
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CBI17951.3 


855 


284 


7 


381 035 


383 836 


2802 


d 


VvWRKYY 


GSVIVTOI 001 286001 


CBI31 897.3 


318 


105 


2 


4 989 461 


4 989 778 


318 


d 


VvWRKY55 


GSVIVTOI 009441 001 


CBI1 9480.3 


960 


319 


18 


8 391 930 


8 393 726 


1797 


d 


VvWRKY23 


GSVIVTOI 022067001 


CBI21 376.3 


843 


280 


7 


16 322 549 


16 324 116 


1568 


d 


VvWRKY36 


GSVIVTOI 029265001 


CBI1 7876.3 


840 


279 


11 


17 821 900 


1 7 823 266 


1367 


d 


VvWRKY12 


GSVIVTOI 0331 88001 


CBI24014.3 


804 


267 


4 


9 363 169 


9 365 026 


1858 


d 


VvWRKY43 


GSVIVTOI 036223001 


CBI35175.3 


915 


304 


14 


8 753 538 


8 756 300 


2763 


e 


VvWRKY32 


GSVIVTOI 021 252001 


CBI30716.3 


837 


278 


10 


3 008 687 


3 010 451 


1765 


e 


VvWRKYS 


GSVIVTOI 01 941 9001 


CBI34393.3 


972 


323 


2 


512 163 


513 533 


1371 


e 


VvWRKY34 


GSVIVTOI 021 765001 


CBI31 090.3 


1266 


421 


10 


10 755 760 


10 759 820 


4061 


e 
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GSVIVTOI 026965001 


CBI40407.3 


1047 


348 


15 


18 957 231 


18 958 817 


1587 


e 


VvWRKY20 


GSVIVTOI 0281 29001 


CBI36956.3 
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7 


4 044 128 


4 045 807 


1680 


e 


VvWRKYS 1 


GSVIVTOI 028823001 


CBI22627.3 


549 


182 


16 


18 360 079 


18 360 711 


633 
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VvWRKY29 


GSVIVTOI 0341 48001 


CBI30536.3 
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8 
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14 830 056 
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Fig. 1. Multiple sequence alignment of the WRKY domain among V. vinifera l/VRKV genes. Red indicates conserved WRKY amino acid domains; green 
indicates zinc-finger motifs; dashes indicate gaps. 'N' and 'C indicate the N-terminal and C-terminal WRKY domain of a specific WRKY gene (this figure 
is available in colour at JXB online). 
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in rice (Wu et al, 2005) and barley (Mangelsen et al, 2008), 
there are no VvWRKY m group-Ill proteins containing a C- 
X7-C-X24-H-X-C zinc-finger motif, perhaps suggesting that 
this is a feature of monocotyledonous species. 

Phylogenetic analysis of V^RKY genes from grape, 
Arabidopsis, and tomato 

To further analyse the evolutionary relationships in the 
Vv WRKY gene family and to help in their classification, a 
total of 210 FP7?^F genes, comprising 70 from Arabidopsis, 
81 from tomato, and 59 from grape, were used to construct a 
phylogenetic tree (Fig. 2). Based on the number of WRKY 
domains and the features of the specific zinc-finger motifs, 
all 59 VvWRKY gQnQS were classified into three main groups, 
with five subgroups in group II (Eulgem et al., 2000). Twelve 
grape WRKY genes {VvWRKY4, VvWRKYll, VvWRKYlS, 
VvWRKYlS, VvWRKY26, VvWRKY28, VvWRKY35, 
VvWRKY39, VvWRKY46, VvWRKY57, VvWRKY58, 
VvWRKY59) with two WRKY domains belong to group I, 
which have a zinc-finger motif of C-X4-C-X22_23-H-Xi-H. 
The other 40 grape WRKY genes with the zinc-finger struc- 
ture of C-X4-C-X22_23-H-Xi-H were assigned to group II, 
which comprised 60% of the total number of Vv WRKY genes. 



The 40 group-II VvWRKY genes are unevenly distributed 
amongst the five subgroups: group Ila (three: VvWRKYQ, 
VvWRKYlO, VvWRKYSO), group lib (eight: VvWRKY2, 
VvWRKYll, VvWRKYSl, VvWRKY38, VvWRKY40, 
VvWRKY45, VvWRKY54, VvWRKY56), group lie 
(15: VvWRKYl, VvWRKY3, VvWRKYS, VvWRKYlS, 
VvWRKYM, VvWRKYM, VvWRKYll, VvWRKY14, 
VvWRKYlS, VvWRKY33, VvWRKY37, VvWRKY44, 
VvWRKY47, VvWRKY49, VvWRKY53), group Ild 
(seven: VvWKYV, VvWRKYll, VvWRKYlQ, VvWRKY13, 
VvWRKY36, VvWRKY43, VvWRKY55), and group He 
(six: VvWRKYS, VvWRKYlO, VvWRKY31, VvWRKY34, 
VvWRKYSO, VvWRKYS 1). In contrast to group I, group-II 
genes have only one WRKY domain. Instead of the C2H2 
pattern, group-Ill genes contain a C2HC zinc-finger motif 
(C-X7-C-X23-H-X1-C) and six of the 59 VvWRKY genes 
{VvWRKY6, VvWRKYlJ, VvWRKY41, VvWRKY41, 
VvWRKY48, VvWRKYSl) belong to this group. Finally, 
Vv WRKYl 7 and Vv WRKY19 do not group within any of the 
other group-II subgroups, possibly because of their appar- 
ently incomplete structures. Detailed information about the 
classification of the genes and the WRKY domains, as well 
as the profile of zinc -finger motifs can be found in Table 1 and 
Supplementary Table S2, respectively. 




Group Mb 

Fig. 2. Phylogenetic tree of genes among grape (red), Arabidopsis (blue), and tomato (black). Filled circle lines are used to cluster the genes with 
similar structures and functions (this figure is available in colour at JXB online). 
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Exon-intron organization of VvWRKY genes 

Insights into the structures of the VvWRKY genes were 
obtained through an analysis of the exon/intron boundaries, 
which are known to play important roles in the evolution of 
multiple gene famihes (Zhang et al. , 2012). As shown in Fig. 3, 
52 of the 59 VvWRKY gsnQS, have two-six exons (seven with 
two exons, 16 with three exons, nine with four exons, nine with 
five exons, and 11 with six exons). The fact that VvWRKY46 
has 17 exons, VvWRKYlS has 11 exons, VvWRKY58 has 
nine exons, VvWRKY38 has eight exons, VvWRKY22 and 
VvWRKYSl have seven exons and VvWRKY? has only one 
exon indicates that both exon loss and gain has occurred dur- 
ing the evolution of the WRKY gene family. This may lead to 
functional diversity of closely related WRKY genes; however, 
this study noted that VvWRKY genes in the same group usu- 
ally have a similar number of exons. The number of exons in 
group I is relatively large, ranging from five to 1 1 ( Fv WRKY46 
with 17 exons was not included because of the possibility of 
special variation), while genes in group II have a relatively 
small number, ranging from one to eight exons, and group III 
from three to five. This exon pattern similarity may be the con- 
sequence of a number of gene duplication events. A further 
analysis indicated that nearly all VvWRKY genes contained 
an intron in their respective WRKY domains (Wu et al., 2005; 
He et al., 2012). The R-type introns are widely exist in major- 
ity of WRKY groups (I, lie, Ild, He, III), the same cases in 
rice and Arabidopsis (Eulgem et al., 2000; Wu et al, 2005). 

Tandem duplication of yy'WRKY genes 

Genomic comparison is a method for rapidly transferring avail- 
able genomic information from a model species to a less-studied 
species (Lyons et al. , 2008; Zhang et al. , 20 12). The current work 
analysed the tandem duplication events of the 59 VvWRKY 
genes on the 19 grape chromosomes (Table 3) according to the 
methods of Holub (2001), where a chromosomal region within 
200 kb containing two or more genes is defined as a tandem 
duplication event. There were 13 VvWRKY genes (VvWRKY9, 
VvWRKYlO, VvWRKYll, VvWRKYlS, VvWRKYM, 
VvWRKYlO, VvWRKYll, VvWRKY24, VvWRKYlS, 
VvWRKY41, VvWRKY41, VvWRKY49, VvWRKYSO) clus- 
tered into six tandem duplication event regions on grape chro- 
mosome 4 (two clusters), 7 (two clusters), 13 (one cluster) and 
1 5 (one cluster) (Table 3). Chromosomes 4 (cluster 1 and cluster 
2) and 7 (cluster 3 and cluster 4) had two clusters respectively, 
indicating a hot spot of WRKY gene distribution. Besides the 
tandem duplication events, 15 segregation duplication events 
were also identified (Fig. 4 and Supplementary Table S3), indi- 
cating that some VvWRKY genes were possibly generated by 
gene duplication. Moreover, the segregation duplication events 
can also provide a reference for the WRKY gene evolutionary 
relationship and functional prediction. 

Synteny analysis of VvWRKY genes 

A substantial number of WRKY genes from the model plant 
Arabidopsis have been systematically investigated (Eulgem 



et al. , 2000; Dong et al. , 2003) and so the current work per- 
formed a synteny analysis of Arabidopsis and grape WRKY 
genes (Fig. 5) to determine whether this might provide some 
functional insights. A total of 66 pairs of syntenic relations 
were identified, including 50 At WRKY genes and 39 Vv WRKY 
genes. Three AtWRKY genes (AtWRKY6, AtWRKY31, 
AtWRKY36) and nine VvWRKY genes (VvWRKY6, 
VvWRKYQ, VvWRKYll, VvWRKYlJ, VvWRKYSl, 
VvWRKYSS, VvWRKYSS, VvWRKY49, VvWRKY54) were 
found to be associated with at least three synteny events 
and interestingly, of these 12 genes, six were in group lib, 
including three AtWRKY genes and three VvWRKY genes 
( VvWRKYSl, VvWRKYSS, VvWRKY54). This may indicate 
a high conservation of group-IIb WRKY genes, which may in 
turn suggest a fundamental function in plant development. 
Detailed results from the comparative analysis are shown in 
Supplementary Table S4. The large number of synteny events 
suggests that many WRKY genes arose before the divergence 
of the Arabidopsis and grape lineages. 

Expression patterns of VvWRKY genes In different 
tissues 

To assess the potential functions of VvWRKY genes during 
grape development, this study investigated the expression 
patterns of all 59 VvWRKY genes in six organs/tissues (roots, 
tendrils, leaves, inflorescences, fruit, and stems). Nearly half 
of the VvWRKY genes showed no significant organ/tis- 
sue related differences in expression, but some clear spatial 
differences were noted (Fig. 6). For example, VvWRKY40 
and VvWRKY45 were expressed at high levels in roots and 
VvWRKY41 and VvWRKYSl expression were particularly 
high in leaves. No Vv WRKY gene was specifically associated 
with inflorescences or fruit. The WRKY genes which showed 
no significant expression difference between tissues are likely 
to play a more ubiquitous role in grape. 

Expression profiles of VvWRKY genes In response to 
drought, salt, and powdery mildew Infection 

The ability of plants to tolerate a variety of abiotic and biotic 
stresses is an essential adaptive feature in changing environ- 
ments. Moreover, the identification and functional analysis 
of genes involved in biotic and abiotic stress signal trans- 
duction pathways is of considerable interest in the context 
of enhancing agricultural productivity. In order to obtain 
insights into the potential roles of all 59 VvWRKY genes in 
stress associated gene expression, this study used 2-year-old 
'Kyoho' {V labrusca x V. vinifera) seedlings growing in pots 
in the greenhouse of Northwest A&F University that were 
exposed to salt and drought stress and monitored expres- 
sion using semiquantitative RT-PCR, 'Shang-24' was used 
for powdery mildew inoculation. A total of four VvWRKY 
genes that could be associated with salt stress {VvWRKY16, 
VvWRKYlS, VvWRKYlS, VvWRKYSS), drought stress 
(VvWRKYS, VvWRKYlS, VvWRKYlS, VvWRKYSS), and 
powdery mildew inoculation {VvWRKYlQ, VvWRKYlV, 
VvWRKY4S, VvWRKYSl), based on semiquantitative 
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Fig. 3. Genomic organization of grape 1/1/RKY genes. (A) Unrooted phylogenetic tree built on tlie basis of 59 complete WRKY-domain proteins in grape; 
details of clusters are shown in different colours. (B) Exon-intron structure of grape 1/VR/<V genes; blue indicates untranslated 5'- and 3'-regions; green 
indicates exons; black indicates introns (this figure is available in colour at JXB online). 
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Table 2. The number of WRKY genes in each phylogenetic group 
from Arabidopsis, rice, grape, and tomato 

VvWRKY17, VvWRKY29, SIWRKY26, SIWRKY27, and SIWRKY49 are 
not included, since tliey could not be placed. 



Gene Phylogenetic group 





1 Ma 




Mb 


lie 


lid 


lie ill 


AtWRKY 


13 4 




7 


18 


7 


9 14 


OsWRKY 


15 4 




8 


15 


7 


11 36 


VvWRKY 


12 3 




8 


15 


7 


6 6 


SIWRKY 


15 5 




8 


16 


6 


17 11 


Table 3. Tandem dupiication events in the 59 VvWRKY genes 


Cluster 


Gene ID 


Chromosome 


Start site 


End site 


number 














1 


VvWRKYQ 


4 






5 247 592 


5 248 886 




VvWRKYlO 


4 






5 265 806 


5 268 041 


2 


VvWRKY12 


4 






9 363 169 


9 365 026 




VvWRKY13 


4 






9 399 944 


9 400 803 




VvWRKYU 


4 






9 409 805 


9 41 1 286 


3 


VvWRKY20 


7 






4 044 1 28 


4 045 807 




VvWRKY21 


7 






4 200 1 60 


4 202 241 


4 


VvWRKY24 


7 






17 794 379 


17 797 240 




VvWRKY25 


7 






17 958 306 


1 7 960 930 


5 


VvWRKY41 


13 






1 716 836 


1 718 836 




VvWRKY42 


13 






1 719 393 


1 720 884 


6 


VvWRKY49 


15 






18 940 954 


1 8 942 1 46 




VvWRKYSO 


15 






18 957 231 


18 957 231 



RT-PCR analysis (Fig. 7A and B and Supplementary Figs 
S1-S3) were selected for further analysis and validation using 
real-time quantitative PCR (Fig. 8A and B). In most cases 
a similar result was seen, including the four genes that were 
induced by drought stress and powdery mildew inoculation 
and three genes (Fvfra^7;6, VvWRKY25, VvWRKY35) by 
salt stress treatment. Only one gene ( Vv WRKY28) showed no 
significant difference in the real-time quantitative PCR. 

Based on the semiquantitative RT-PCR data (Fig. 7A), 
VvWRKY genes tend to be downregulated to a greater 
degree by salinity stress than by drought stress. VvWRKYll, 
VvWRKYM, VvWRKYlS, VvWRKY26, VvWRKY28, 
VvWRKYSl, VvWRKY32, VvWRKY39, VvWRKY46, and 
VvWRKY48, which all showed clear downregulation by the 
salt stress treatment, were upregulated to varying degrees 
by the drought stress treatment, indicating distinctly dif- 
ferent regulatory networks existence. The reaction time for 
the expression pattern changes was another focus of this 
study and some genes, such as VvWRKYl and VvWRKYSl, 
showed altered expression at an early time point (Ih under 
salt stress, 24 h under drought stress), while others, includ- 
ing VvWRKY57 and VvWRKY59 showed upregulated 
expression at a relatively late time (48 h under salt stress, 
144 or 168 h under drought stress). On the other hand, 
some genes showed an early downregulation but subsequent 



upregulation (VvWRKY2, VvWRKY27, VvWRKY49): for 
example, Vv WRKY2 was downregulated at 1 h after the onset 
of the salt stress treatment, but was upregulated after 48 h for 
example. 

The association of the VvWRKY gQne?, with biotic stress 
was investigated using infection with powdery mildew, a 
severe disease worldwide (Zhang et al., 2012). The semi- 
quantitative RT-PCR based expression profiles (Fig. 7B) 
of VvWRKY4, VvWRKYl 1, VvWRKY12, VvWRKY23, 
VvWRKY24, VvWRKY30, VvWRKY31, VvWRKY32, 
VvWRKY33, VvWRKY38, VvWRKY40, VvWRKY42, 
VvWRKY43, VvWRKY46, VvWRKY47, and VvWRKYSl 
indicated a potentially negative effect on grape powdery mil- 
dew resistance, while VvWRKY?, VvWRKY8, VvWRKYQ, 
VvWRKYl 3, VvWRKYM, VvWRKY25, VvWRKY26, 
VvWRKY27, VvWRKY34, VvWRKY35, VvWRKY36, 
VvWRKY37, VvWRKY39, VvWRKY41, VvWRKY44, 
VvWRKY45, VvWRKY48, VvWRKY49, VvWRKYSO, 
VvWRKY52, VvWRKY55, and VvWRKY58 showed the 
opposite pattern. VvWRKY12, VvWRKY23, VvWRKY40 
VvWRKY43, and VvWRKY46 showed significant down- 
regulation during early infection (6-72 h or to 96 h), while 
VvWRKY34, VvWRKY35, and VvWRKY36 were upregu- 
lated during a later stage (24-120 h), possibly indicating a 
relationship between expression pattern and responding time. 

Expression profiles of VvWRKY genes to liormone 
treatments 

Plant hormones such as ABA, SA, MeJA, and ethylene have 
well-established roles in modulating plant signalling networks 
(Fujita et al. , 2006). In this study, hormone treatments resulted 
in a wide variety of VvWRKY gene expression profiles (Fig. 7C 
and Supplementary Figs S4-S7). A total of 31 VvWRKY genes 
showed different degrees of downregulation by the ABA treat- 
ment while 14 were upregulated. Similarly, 19 VvWRKY genes 
were downregulated and 18 were upregulated expression fol- 
lowing SA treatment. However, the expression profiles result- 
ing from MeJA and Eth treatments were distinct from those 
modulated by ABA and SA and a substantially greater number 
of upregulated genes were observed: 37 were upregulated by 
MeJA and 11 were downregulated, while 37 were upregulated 
and six were downregulated by Eth treatment. These expres- 
sion variations indicate that the Vv WRKY gene family possible 
is collectively regulated by a broad set of hormonal signals. 

Discussion 

Many WRKYfamily genes play important roles in diverse plant 
developmental and physiological processes (Du and Chen, 
2000; Eulgem et al, 2000), including embryogenesis (Lagace 
and Matton, 2004), seed coat and trichome development 
(Johnson et al. , 2002), and leaf senescence (Miao and Zhetgraf, 
2007; Zhou et al, 2011), as well as various plant abiotic and 
biotic stress responses. This current study describes the iden- 
tification of 59 VvWRKY genes from grape, together with an 
analysis of their structure, evolutionary history, and expression 
pattern diversity with respect to biotic and abiotic stresses. 
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Fig. 4. Chromosome distribution and synteny analysis of grape WRKY genes. Chromosomes 1-19 are shown with different coiours and in a circuiar 
form. The approximate distribution of each VvWRKY gene is marked with a short red line on the circle. Coloured curves denote the details of syntenic 
regions between grape WRKY genes (this figure is available in colour at JXB online). 



Identification and annotation of VvWRKY genes 

The publicly available collection of V vinifera ESTs were 
used to confirm the identity of the 59 VvWRKY gwQS that 
were initially identified amongst the 30 434 annotated grape 
genes using BLASTP. Of the 59 VvWRKY gQnQ?,, seven were 
without supporting EST data, indicating that they may not 
be expressed during grape development. However, semi- 
quantitative RT-PCR analysis was used to confirm that all 
59 VvWRKY gQnQ?,, including the seven genes without EST 
support, were indeed expressed and had putative functions 
in various aspects of grape biology, indicating all these 59 
VvWRKY gQnQS, were putative WRKY gQrvQi, in grape. 

Structural conservation and divergence of 
VvWRKY genes 

Comparative genomic analysis is an effective method for study- 
ing gene structures and so this study assessed the conserved 
structural domains of the predicted grape WRKY proteins. 
Multiple sequence alignments revealed that five VvWRKY pxo- 
teins (VvWRKYS, VvWRKYB, VvWRKYM, VvWRKYlJ, 
VvWRKY24) had sequence variations in their WRKY 
domain and three {VvWRKY7, VvWRKYl?, VvWRKY46) 



had variations in their zinc-finger motif In previous studies of 
Arabidopsis transcription factors, it was found that 

the binding-site preferences of the WRKYGQK motif depend 
on the DNA sequences adjacent to the TTGACY core motif 
(Ciolkowski et al., 2008). In the genes recognized by the pro- 
teins with a variation to the WRKY domain, the TTGACY 
core motif may exhibit an altered structure and function, and 
so WRKY genes without the WRKYGQK motif may recog- 
nize binding sequences other than the W-box element ((C/T) 
TGAC(C/T)). In tobacco, the NtWRKY12 protein, which 
contains a WRKYGKK motif, recognizes the downstream 
binding sequence TTTTCCAC, which is substantially differ- 
ent from the W-box (van Verk et al. , 2008). Moreover, the soy- 
bean WRKY proteins GmWRKY6 and GmWRKYll, which 
have a WRKYGKK motif, do not bind normally to the W-box 
(Zhou et al, 2008). It therefore seems that variations in the 
WRKYGQK motif influence the normal interaction of WRKY 
genes with downstream target genes, and so it would be inter- 
esting to investigate the functions and binding specificities of 
VvWRKYS, VvWRKYB, VvWRKYM, VvWRKYl?, and 
VvWRKY24. Moreover, as far as is known, nothing is known 
about the effect of zinc-finger motif changes. It is possible that 
the variations of the WRKY domain and zinc-finger motif could 
influence the classification of the VvWRKY genes reported here. 
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Fig. 5. Synteny analysis of grape and Arabidopsis WRKY genes. The chromosomes of grape and Arabidopsis are depicted as a circie. The approximate 
distribution of each AtWRKY gene and VvWRKY gene is marl<ed with a short red line on the circie. Coioured curves denote the details of syntonic regions 
between grape and Arabidopsis WRKY genes (this figure is available in colour at JXB online). 
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Fig. 6. Expression pattern of 59 VvWRKY genes in organs/tissues of 'Kyoho' {V. labrusca x V. vinifera); Actin 1 (GenBank accession number AY680701 ) 
was used as an internal control. Lanes: 1 : root, 2: tendril, 3: leaf, 4: inflorescence, 5: fruit, 6: stem. The experiments were repeated three times and the 
results were consistent. 
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One example of this is VvWRKYlJ, which has a divergent zinc- 
finger sequence, leading to a nebulous classification in the phy- 
logenetic tree. It remains to be determined whether the observed 
sequenced variations in the conserved domains affect the func- 
tion or the expression patterns of the regulated gene targets. 

Exon-intron structural diversification also plays an impor- 
tant role in the evolution of many gene families and exon- 
intron gain or loss may be caused by the rearrangement and 
fusions of different chromosome fragments (Xu et a!., 2012; 
Guo et a!., 2013). The current study provides an example of 
such diversification in the form of a WRKYgmc ( VvWRKYT) 
with only one exon, while other genes in the same phyloge- 
netic group (group lib) have four or five exons. Moreover, 
VvWRKY46 has as many as 17 exons, while VvWRKY4, 
which is similar to VvWRKYl?, only has six exons. 

The evolutionary relationship of VvVMRKY genes 

The size of the grape WRKY gmc family (59) is small com- 
pared to that of other experimental model plants such 
as Arabidopsis (72) and rice (96). The WRKY genes from 
grape, Arabidopsis, rice, and tomato align within distinct 



phylogenetic groups (Table 2) and it is apparent that vari- 
ations in the number of WRKY genes in group III are the 
primary cause of the diversity of WRKY gene family size. 
Interestingly, previous studies described group-Ill WRKY 
genes as being a newly defined group, as well as being the 
most dynamic group with respect to gene family evolution 
(Zhang and Wang, 2005). Therefore, a key role of group-Ill 
WRKY gQviQ^ in plant evolution may exist. 

Gene duplication events play a major role in genomic 
rearrangements and expansions (Vision et al. , 2000) and are 
defined as either tandem duplications, with two or more genes 
located on the same chromosome, or segmental duplications, 
with duplicated genes present on different chromosomes (Liu 
et al. , 2011). The large number of gene duplication events for 
grape (Figs 4 and 5) will help aid future gene function pre- 
diction and evolution analysis. Whole-genome duplication 
events (y, (3, a) are a common phenomenon in angiosperms 
(Zhang and Wang, 2005) and often lead to gene family expan- 
sion (Cannon et al, 2004). Ling et al (2011) reported that in 
cucumber {Cucumis sativus) CsWRKY family, a divergence 
generated in the number of group-Ill I^i^^F genes resulted 
from different style of duplication events that occurred after 
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the divergence of the eurosids' groups I and II (110 Mya). which may be caused by a different pattern of duplication 
For some species in the eurosids' group I (cucumber, soybean, events. This idea is consistent with the current results concem- 
and grape), the number of group-Ill WRKY genes is small, ing the group-Ill VvWRKY genes. The group-Ill AtWRKY 
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genes with normal tandem duplication events (AtWRKY63, 
AtWRKY64, AtWRKY66, AtWRKY67) show evidence of 
large-scale duplication. According to Zhang (2003) not all 
duplication events are stable, but instead can be fixed or lost in 
the population due to selection pressure and evolution. If the 
duplication events happens at a favourable time and in genes 
that are highly expressed, they will most likely be retained, but 
other duplication events that are not useful to the organism 
because of functional redundancy, or even a negative effect 
on plant development, will either be deleted from the genome 
or become very will diverge. The group-Ill WRKY gtrve diver- 
gence may reflect such a gene evolutionary selection. 

VvWRKY genes function in abiotic and biotic stresses 

There is considerable evidence that WRKY genes play crucial 
roles in responses to abiotic and biotic stress-induced defence 
signalling pathways (Qiu et ciL, 2004; Chen et al., 2012; 
Ishihama and Yoshioka, 2012). From an applied perspective, 
the identification of WRKY genes with potential value in stress 
resistance improvement of grape would likely benefit from tar- 
geting such genes that are part of abiotic and biotic stress- 
response networks. The semiquantitative RT-PCR expression 
profiles generated in this study (Fig. 7) revealed different 
expression patterns (upregulation and downregulation) for 
each VvWRKY gene under specific treatments, thus provid- 
ing a useful resource for future gene expression and functional 
analyses. Dong and coworkers (2003) showed that nearly 70% 
of the At WRKY genes are differentially expressed in response 
to microbial infection or SA treatment. Consistent with these 
previous studies, the current results show that 70-90% of 
VvWRKY genes are differentially expressed following various 
abiotic and biotic stress treatments, highlighting the extensive 
involvement of WRKY genes in environmental adaptation. 

The roles of many Arahidopsis WRKY genes in plant abi- 
otic stress responses have been extensively studied recently. 
For example, AtWRKY25 and AtWRKY33 regulate plants 
adaptation to salinity stress through an interaction with 
their upstream or downstream target genes (Jiang and 
Deyholos, 2009). VvWRKY26 shares 76% sequence similar- 
ity with AtWRKY25 and 71% similarity with AtWRKY33, 
but in the current study, the change in expression profile 
of VvWRKY26 as a consequence of the various treat- 
ments appHed was not apparent. It is possible that although 
VvWRKY26, AtWRKY25, and AtWRKY33 are segregation 
duplication genes, their function may differ in different plant 
species. The presence or absence of the regulatory elements in 
the duplicated genes can have an important consequence on 
subsequent divergence of gene function for an explanation. 

VvWRKYl (Marchive et al, 2007) and VvWRKY2 (Mzid 
et al. , 2007) (named in this work Vv WRKY53 and Vv WRKY4, 
respectively) have also been isolated and functionally analysed. 
VvWRKY2 is known to activate the promoter of VvC4H, 
which is involved in the lignin biosynthesis pathway and cell- 
wall formation (Guillaumie et al. , 2010), and it hkely plays an 
important role in resistance or tolerance to biotic and abiotic 
stress in plants. However, VvWRKY4 appears not to respond 
to powdery mildew inoculation, suggesting that it may involve 



in resistance other abiotic and biotic stresses other than 
powdery mildew. On contrary, VvWRKY53 is significantly 
upregulated at 24 h post inoculation, which appears to share 
similar inoculation response with Vv WRKYl as reported ear- 
lier (Marchive et al, 2007), thus suggesting that VvWRKY53 
may play a role in eliciting resistance response during early 
stage of infection. 

The temporal and spatial diversification of WRKY 
gene expression is widespread, which is important for 
gene function analysis. The VvWRKY genes VvWRKY40, 
VvWRKY45, VvWRKY42, and VvWRKY52 had a higher 
expression in certain grape organs/tissues suggesting diver- 
gent roles for these genes in grape development. With 
regard to temporal diversification, this work considered the 
response time to the different treatments since previous stud- 
ies had shown that some WRKY genes respond to drought 
stress at an early stage, such as Gs WRKYl 8, which peaked 
at 0.5 h after drought stress treatment (Ling et al., 2011). 
The current data show the opposite pattern, since about 
half of the VvWRKY genes {VvWRKY3 and VvWRKY35) 
selected for real-time quantitative PCR showed a peak of 
expression at 144 h after treatment and one {VvWRKY28) 
peaked at 168 h after drought stress treatment, indicating 
that a certain response time is needed for VvWRKY genes to 
respond to drought stress. The upregulation of the expres- 
sion of some WRKY genes in response to powdery mildew 
inoculation showed a similar delay in response time, but 
shorter than that caused by the drought stress treatment 
(12 h post treatment). 

WRKY genes that are components of plant biotic stress regu- 
latory networks have a complex response pattern. For example, 
Arahidopsis AtWRKY33 (Lippok et al , 2007) and AtWRKYlS 
(Chen and Chen, 2002) can positively modulate defence-related 
gene expression and improve disease resistance, while some 
negative regulatory elements may prevent the overexpression 
of AtWRKY33 and AtWRKYlS and are detrimental to plant 
growth. Moreover, other WRKY genes, such diS At WRKY7{¥Jim 
et al, 2006) and AtWRKY48 (Xing et al, 2008) have an imme- 
diate negative effect in the plant defence response. In these cur- 
rent studies, although VvWRKY4, VvWRKYll, VvWRKY12, 
VvWRKY23, VvWRKY24, VvWRKY30, VvWRKY31, 
VvWRKY32, VvWRKY33, VvWRKY38, VvWRKY40, 
VvWRKY42, VvWRKY43, VvWRKY46, VvWRKY47, and 
VvWRKYSl may have a negative effect on grape powdery 
mildew resistance and VvWRKY?, VvWRKYS, VvWRKYQ, 
VvWRKYl 3, VvWRKYl 4, VvWRKY25, VvWRKY26, 
VvWRKY27, VvWRKY34, VvWRKY35, VvWRKY36, 
VvWRKY37, VvWRKY39, VvWRKY41, VvWRKY44, 
VvWRKY45, VvWRKY48, VvWRKY49, VvWRKYSO, 
VvWRKY52, VvWRKY55, and VvWRKY58 may have a posi- 
tive effect, the possible interactions between two or more genes 
and regulatory mechanisms remain to be resolved. 

Phytohormones are critically important in coordinating reg- 
ulatory networks and the signal transduction pathways associ- 
ated with external cues. ABA is known to function in signalling 
in some stressful environments (Huang et al, 2008), and SA, 
JA, and Eth play important roles in biotic stresses (Fujita 
et al. , 2006), while MeJA responds to some biotic stresses and 
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wounding (Huang et al., 2008). The current results show that 
the majority of VvWRKY genes, showed significantly upregu- 
lated expression following treatment with MeJA and Eth, while 
ABA and SA treatments had the opposite effect. It has been 
reported that AtWRKY33 can act as a positive regulator in 
Eth-mediated defence signalling against necrotrophic patho- 
gens and as a negative regulator in SA-mediated responses 
for some biotrophic bacterial pathogens (Zheng et al., 2006), 
which is the same trend seen in this study. 

In conclusion, WRKY gene expression is influenced by a 
broad range of abiotic and biotic stress resistances and also 
responds to hormonal signals, indicating that they play a 
key role in signal transduction in plant resistance regulation. 
The role of WRKY genes on abiotic and biotic stress regula- 
tory networks is systematic and complex, both between two 
WRKY proteins (Xu et al. , 2006) and their downstream or 
upstream targets (Jiang and Deyholos, 2009). 
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